Geometrically-controlled twist transitions in nematic cells 
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We study geometrically-controlled twist transitions of a nematic confined between a sinusoidal 
grating and a flat substrate. In these cells the transition to the twisted state is driven by surface 
effects. We have identified the mechanisms responsible for the transition analytically and used exact 
numerical calculations to study the range of surface parameters where the twist instability occurs. 
Close to these values the cell operates under minimal external fields or temperature variations. 



PACS numbers: 61.30.Cz, 61.30.Gd 

Liquid-crystal (LC) displays consist of a LC layer con- 
fined between two surfaces that impose a preferred orien- 
tation of the average molecular direction n(r). Current 
nematic displays rely on voltage-induced reorientation of 
the director within the bulk LC layer. Recently, textured 
surfaces (on scales of the order of 1/im) were designed for 
patterned alignment of LC Jl]-||] , opening possibilities for 
an improved performance of LC cells. In view of the very 
rich behavior of LCs even on homogeneous substrates jj| 
and of simple fluids on structured surfaces such tech- 
nological advances beyond a trial and error procedure re- 
quire theoretical guidance. As a first step in this direction 
we establish the minimal model appropriate to the theo- 
retical description of patterned LC cells and develop an 
efficient method of solving it. This study contributes to 
the broader perspective of how well-defined structures on 
solid surfaces can be imprinted on adjacent soft matter. 
The particular softness of LCs leads to the expectation 
of very pronounced effects. 

The microscopic description of LC surfaces and inter- 
faces is complex. The degree of nematic order, biaxility, 
etc. vary in the interfacial region over molecular distances 
while the director may vary over macroscopic distances. 
In LC displays, the characteristic distance over which 
the director varies is set by the cell dimensions or by the 
electric correlation length Q. Under most experimen- 
tal conditions the length over which the director varies 
is of the order of /jm and a macroscopic or elastic the- 
oretical description is adequate The macroscopic 
(second order) elastic free energy of bulk nematics was 
established more than 40 years ago |J but the status of 
the surface contributions in the weak anchoring regime 
is still controversial 1 10 Tl|] . 

In the following we propose and minimize an elastic 
free energy for patterned LC displays that includes bulk 
and surface terms. We consider the twist cell proposed 
in a recent experiment, where the nematic is confined be- 
tween a flat and a sinusoidal grating surface pj. In the 
experiment a voltage-controlled twist (VCT) effect highly 



sensitive to the surface properties of the cell (grating ge- 
ometry and anchoring strength) has been reported, for 
gratings on the scale of tenths of fim. Berreman and 
de Gennes [Q first considered grating surfaces to explain 
azimuthal anchoring by elastic effects only. Faetti Jl3"| , 
and more recently Fournier and Galatola [03, general- 
ized the effective azimuthal anchoring energy by intro- 
ducing local anchoring at the grating surface. Barbero 
and Durand [ fl5| also considered grating surfaces with 
characteristic lengths comparable to the nematic correla- 
tion length. They used the Landau-de Gennes free energy 
to describe the induced quasi-melting caused by the 
rough surface. Over the last few years there has been a 
considerable surge of interest in the influence of grating- 
like surfaces on the structural properties of LCs because 
new techniques enabled the manufacture of controlled 
undulated surfaces, allowing for meaningful comparisons 
between theory and experiments. Among these new sys- 
tems the aforementioned twist displays exhibit excellent 
viewing angle characteristics [Q, that are important for 
technological applications. 




FIG. 1. Nematic director field (projected onto the xz plane) 
for increasing anchoring strengths W\. (a) Distorted director 
field, for strongly negative anchoring strength W\. (b) Upon 
increasing Wi a first instability leads to a nearly uniform di- 
rector field, (c) For larger Wi the director field follows the 
orientation of the surface, (d) A second instability induces a 
twisted director field. 
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In spite of its importance, a theoretical study of twist 
cells covering the whole range of operating parameters is 
still lacking. In order to explore how variations in the 
geometry of the cell can induce the reorientation of the 
bulk nematic director we consider the total free energy 
of the LC cell as the sum of a bulk elastic free energy Fb 
and a surface contribution F s . The former is the Frank 
elastic free energy |?J 

Fh = \ J v {Kll(y ' + K22[U ' (V X n)]2 

+ #33 [n x (V x n)] 2 }d 3 r (1) 

where K\\, K22, and -K33 are the elastic constants associ- 
ated with splay, twist, and bend distortions, respectively. 
F s includes the anchoring energy for which we adopt the 
Rapini-Papoular form pi 

Fs = ^j (n.„)Vr+^»/ (n.,)^r. (2) 

J Si <J S2 

These integrals run over the two cell surfaces and Wi, 
i = 1,2, is the corresponding anchoring strength that 
characterizes each surface, v is the local unit vector nor- 
mal to the surface. For negative Wi this energy contri- 
bution favors normal surface orientation of the nematic 
while positive Wi favor planar (degenerate) orientation at 
the surface. For inhomogeneous substrates in the weak 
anchoring regime one may have to include a surface elas- 
tic term associated with the saddle-splay distortion jl0| . 
Its contribution will be considered later. 
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FIG. 2. The twist instability (numerical results). The twist 
order parameter N = D~ 3 f^, njjd 3 r as a function of the 
reduced anchoring strength Wi = WiD/Kn for different 
groove depth A/D. To favor twist configurations we used 
Kn = K33 — 2X22- N vanishes linearly at Wi c . 

The total free energy, F t — Fb + F s , is a functional of 
the two angles a(r) and 0(r), characterizing the nematic 
director n = (cos a cos 0, cos a sin 0, sin a) . The cells con- 
sidered in this letter consist of a nematic LC confined 
between a sinusoidal grating {z = D + Asinqx, where 
q = 2n/L and A is the groove depth) and a flat sub- 
strate (z = 0) (see Fig. 1). In Ref. jl| the surfaces were 



treated so that the grating induces (weak) homeotropic, 
i.e., normal anchoring with respect to the local surface 
direction, while the flat surface induces (strong) homoge- 
neous uniaxial anchoring in the direction perpendicular 
to the grooves. In the following calculations we keep the 
homogeneous strong anchoring condition at the flat sur- 
face (a(x, z = 0) = <f>(x, z = 0) = 0) but we vary the 
anchoring strength W\ (including its sign) at the grating 
surface. Finally, we use periodic boundary conditions 
a{x = 0, z) = a(x = L, z) and <p{x = 0, z) = <p(x = L, z). 

Figure 1 illustrates the behavior of a confined LC with 
a small twist elastic constant as the anchoring strength 
W\ increases. For large negative values of Wi, the grat- 
ing surface induces normal orientation of the nematic. 
Upon increasing W\ a first instability occurs so that the 
director field becomes nearly uniform resulting from the 
competition of two effects: the anchoring energy favor- 
ing homeotropic anchoring, and the elastic energy favor- 
ing homogeneous alignment at the grated surface. The 
critical anchoring strength for a rectangular cell D x S, 
where a(x, z) = az/D, is found easily. For small a, the 
total free energy is F t = (S/2)(Wi + Kn/D)a 2 , where 
S is the area of the flat cell surface. The director field 
is uniform when the coefficient of a 2 is positive, i.e., for 
W\ > —Kn/D. In the limit of small groove depth A 
the corrections arising from the grating are found by 
first-order perturbation analysis about the rectangular 
cell. When W\ becomes positive, the nematic orienta- 
tion changes continuously to follow the sinusoidal shape 
of the boundary. Beyond a certain threshold, the bulk 
bending energy of this deformation is comparable with a 
bulk twist deformation and a new instability occurs. Ow- 
ing to the degeneracy of the planar anchoring, the sinu- 
soidal and the twisted configuration have approximately 
the same surface free energy. 
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FIG. 3. The onset of the twist instability occurs for an- 
choring strengths Wic which depend on the cell geometry 
(A, q = 2-k/L). Here the contour lines for W\ c = 0.5, 1, 5, 10, 
and 15 are shown for Kn = K33 — 2X22- These results 
are obtained analytically based on Eq.(^). The vertical and 
horizontal line indicates systems with L — D or A = -D/10, 
respectively (c.f., Fig. 4 (b),(c)). 
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The twist instability depends strongly on the geome- 
try of the cell. In Fig. 2 we plot the numerical results for 
the twist order parameter N — D~ 3 J y n 2 d 3 r as a func- 
tion of the anchoring strength W\ for cells with different 
groove depth A. 

The free energy of the sinusoidal configuration F S i„ is 
estimated by using the variational ansatz 



. sinh Xz 

a s in(x,Z) = BAq cos qx — — — 
sinh XD 



(3) 



This function respects the boundary conditions and has 
two parameters B and A. If the nematic director field is 
nearly uniform the dimensionless parameter B ~ 0. On 
the contrary, if the director is parallel to the sinusoidal 
surface, B = 1. A -1 defines the distance over which the 
nematic deformation is influenced by the grating surface. 
For small groove depth A/D the total energy is 



K 



S 



.„ „ -(Aq) 2 [Y(X)B 2 + - Bf] 



(4) 



where Y{\) = Ki 1 X 2 F 1 (X) + K 33 q 2 F 2 (X), with 
Fi = (cosh Xz/ sinh \D) 2 dz and F 2 = 

§q (sinh Xz / sinh XV) 2 dz. Minimization with respect to 
the variational parameters yields A = q yj K33/K11 and 
B = W\/(W\ + Y). The total free energy follows as 



S 



YWi 



(•5) 



F. 



(c) _ 
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22 



2D 



(9) 



In Fig. 3 we plot the results obtained from solving this 
equation for different cell geometries (A,q = 2w/L). To 
favor twist configurations we used K\\ = K 33 = 2K 22 . 

In order to overcome the limitations of the above an- 
alytic analysis extensive numerical calculations of cells 
with a wide range of surface parameters (groove depth, 
pitch, and anchoring strength) have been carried out. 
The calculation of the twist transition requires a numeri- 
cal procedure capable of evaluating the bulk and surface 
contributions to the free energy very accurately. Ow- 
ing to the geometrical pattern of the grating surface this 
turns out to be a challenging numerical problem. We 
used finite-element discretizations of the functions a(x,z) 
and <f>{x, z) (dividing the space into small triangles where 
the functions are approximated linearly) and found the 
minimum of the free energy by standard numerical tech- 
niques. Non-uniform adaptive meshes have been used: 
finer meshes were required close to the grating surface 
where the fields vary more rapidly Jl7[ |. In addition, a 
finite-element triangulation that approximates as closely 
as possible the geometrical boundary of the cell was re- 
quired since the critical anchoring strength depends sen- 
sitively on these surface terms. 



For Wi/Y ^> 1 and XD > 1 we recover the Berreman 
energy 



-(Aq) 2 qy/K 33 K n 



(6) 



which is independent of the cell size D. However, when 
the flat surface approaches the grating one finite size 
effects come into play. In this case, for AD <C 1, 
Y w [K 11 +K 33 {qD) 2 /3]/D. 

In order to estimate the critical anchoring strength 
Wic, i.e., the threshold for the twist instability we take 
4>{x, z) = <f>z/D and use an expansion for the total energy 
in powers of <f>: 



F t [a, 



F, 



(2) 



(7) w. 



The sinusoidal deformation a s 



0) minimizes 

the free energy if the coefficient F^ [a S m] is positive. 
Supposing that K 22 /K\\ ~ (A/D) 2 is small one has 



(with ■ 



F 



(2) 



_ S -(Aq) 2 ^ 

2D A y H ' Y 



Wi 



(8) 



The twist instability occurs for the threshold value given 
by _F t [a s in] — 0, or equivalently when the energy of the 
sinusoidal deformation equals 




qD 

FIG. 4. Critical anchoring strength W\ c (A,q) for 
Ku — K33 — 2K22 as in Fig. 3 obtained from an exact nu- 
merical solution, (a) Contour lines W\ c = 0.5,1,5,10, and 
15. (b)W la (A,L = D). (c)W lc (A = 0AD,q). 
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These numerical results are shown in Fig. 4. From this 
figure one infers that W\ c diverges for a set of the geomet- 
rical parameters of the cell. This set defines a range of 
parameters for which the nematic orientation will never 
twist. Outside this region, for larger values of A and q, 
a twisted nematic director may occur for small values of 
W lc = W lc D/K n . 

For completeness, we consider the elastic term asso- 
ciated with the saddle-splay distortion neglected in the 
previous analysis. This term is usually written in the 
form 

F ss = K 2A [ V • [(n • V)n - n(V • n)] d 3 r (10) 
Jv 

which reveals its surface nature. The value of the surface 
elastic constant K 2 4 is bounded for stability reasons. In 
fact, the bulk free energy in Eq.(jl]) is the sum of quadratic 
terms and it is well defined only if the bulk elastic con- 
stants K\x, K22, and K33 are positive. A stability analy- 
sis of the surface energy in Eq. ( |l0| ) yields the constraint 
< if 24 < min(ifii, K22) (jj]. In Fig. 5 we plot the nu- 
merical solutions for a (x, z) for physical and unphysical 
values of if 24- By varying if 24 (between and if 22) we 
found that the thresholds W\ c for the twist instability 
are not significantly affected by this term |[9| within the 
physically relevant range. 




FIG. 5. Numerical solutions for the system considered 
above (with A/D = 0.1, qD = 2tt, Wi = 2, and 
Kn — K33 — 2X22)- The function \a(x, z)\ is represented 
by a color code in which blue and red correspond to a = 
and a = max|a|, respectively, for physical (left: K24, = 0) 
and unphysical values of -K24 (right: K24 = 2Kn). The un- 
physical solution on the right exhibits strong, spurious defor- 
mations close to the boundary leading to a divergent negative 
free energy that is bounded numerically due to the mesh dis- 
cretization. The left panel corresponds to the bent sinusoidal 
configuration shown in Fig. 1(c). 

Beyond the obvious technological importance of the 
twist cell - it has been shown to possess an electro-optic 
response far less dependent on viewing angle than other 
LC display configurations - the system turns out to 
be very interesting also from a theoretical point of view: 
for a given anchoring strength, the twist transition is 



driven in the absence of an external field by the surface 
morphology (A/D,qD). Close to the surface transition, 
an arbitrarily small external field or temperature varia- 
tion will be sufficient to induce the reorientation of the 
nematic. 

An obvious extension of this work includes application 
of a voltage between conducting substrates. This prob- 
lem, however, requires special care. At constant voltage 
the system is no longer isolated and the minimum prin- 
ciple for the total free energy (with electric and elastic 
terms) does not apply. The free energy is minimal with 
respect to the nematic director field and maximal with 
respect to the elecric potential ^0|. A generalization of 
the numerical method described above may not converge 
and a numerical solution of the Lagrangian differential 
equations appears to be required. 

PP acknowledges the support of the Fundagao para 
a Ciencia e Tecnologia (FCT) through Grant No. 
SFRH/BPD/5664/2001. 
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